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. Abstract: Assuming (i) the seesaw mechanism works for neutrino masses and mixings, {ii) dark 
I matter is absolutely stable due to unbroken U{l)x dark gauge symmetry, and [iii) the singlet fields 
{H^H and Nr) are portal to the dark sector, we construct a simple model which is consistent with all the 
\ cosmological observations as well as terrestrial experiments available as of now, including leptogenesis, 
extra dark radiation of ~ 6% (resulting in A^eff = 3.106 the effective number of neutrino species), Higgs 
inflation, small and large scale structure formation, and current relic density of scalar DM {X). The 
' electroweak (EW) vacuum of this model is stable up to Planck scale for mn = 125 GeV without any 
. other new physics. The Higgs signal strength is equal to one as in the SM for unbroken U{\)x case 
: with a scalar dark matter, but it could be less than one independent of decay channels if the dark 
ff^ , matter is a dark sector fermion or if U(l)x is spontaneously broken, because of a mixing with a new 
neutral scalar boson in the models. Detailed study of Higgs properties at the LHC would shed light 
on the models described in this work. 
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1 Introduction 

The SM has passed many tests from various experiments from atomic physics scale up to a couple of 
TeV energy range. Still it is well known that the SM has to be extended in order to accommodate the 
neutrino masses and mixings, baryon number asymmetry (BAU) and cold dark matter (CDM) of the 
universe. The most economic explanation for the 1st two problem would be leptogenesis [1, 2], whereas 
there are many models for CDM in particle physics [3-6]. 

For CDM physics, one of the puzzles is how CDM can be absolutely stable or very long lived. 
If unstable, the lifetime of CDM should be far longer than the age of the universe, say r > 10^^"^*^ 
sec [7]. Otherwise its decay would produce too much A(7)-ray or neutrino flux to match observation. 
Still this lower bound of the CDM lifetime is far less than the lower bound on the proton lifetime, the 
reason of which still remains one of the mysteries in particle physics. 

The required longevity of the DM lifetime can be guaranteed by a symmetry. If the symmetry is 
global, it can be broken by gravitational effects, and there can be dangerous operators suppressed by 
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Planck scale (Mp), such as 



■-decay 



^ <\)xF^uF^'' for bosonic DM 



igf- {Ipf-Li) for fermionic DM ^ 



(1.1) 



where \^ = O(e^) and ~ C'(l) are the couplings associated with the operators, (j)x and ■i/'x are 
bosonic and fermionic dark mater candidates, and In and H are the SM lepton and Higgs, respectively. 
In this case, dark matter can not be stable enough unless < 0(10) keV for bosonic CDM (/>, or 
< 0{1) GeV for fermionic CDM ip. It may be possible to have such a light dark matter though 
it may not be theoretically so natural. Axion or keV scale sterile neutrinos are good examples of DM 
whose longevity is guaranteed by some global symmetries. On the other hand, the argument above 
implies that it is highly unlikely that an EW scale CDM is long lived or stable due to some global 
symmetries. 

Contrary to global symmetry, local symmetry other than the SM gauge group often appears in 
theories beyond the SM, and would guarantee the absolute stability of dark matter if it is unbroken ^. 
For example, gauge groups in superstring theory have very large ranks, e.g. SO(32) or x E'^. At low 
energy, these gauge groups may be broken into a product of the SM gauge group {SU{'i)c x SU{2)l x 
U{1)y) and some other group, the latter of which may be able to play a role of dark gauge group 
we consider in this paper. The presence of a unbroken extra local symmetry implies the existence of 
massless gauge boson(s). Since it is a carrier of a long range force ^, the massless dark gauge boson(s) 
could have significant effects on structure formation via self-interactions of dark matter [16-18]. On 
one hand, the dark gauge interaction is highly constrained by various properties of small and large 
scale dark matter halos [19]. On the other hand, it can provide a solution to small scale puzzles of 
the collisionless CDM scenario (e.g. cored density profiles [20-23] and low concentrations of massive 
sub-halos [24]) without conflicting with constraints from large scale structure [25]. The massless dark 
gauge boson(s) could also contribute to the radiation density of the Universe in addition to thermal 
relic neutrinos of 3 species. Recent WMAP 9-year data analysis showed that the number of relativistic 
degrees of freedom is [26] 

N°^' = 3.84 ± 0.40 at 68% CL (1.2) 

Although it is consistent with the case of three active standard model neutrinos only (A'^f^ = 3.04), 
some amount of extra radiation is still allowed and it could be from either light sterile neutrino [27] or 
hidden photon [28] or axion [29-36]. There are considerable amount of literatures on these possibilities. 

Meanwhile, dark sector can communicate with the SM sector via Higgs portal interactions {H^H) 
which is quite often used in the dark matter physics [37-42] (see also [43] where DM produced from SM 
particles via kinetic mixing and Higgs portals was analyzed). Another possible portal interaction can be 
provided by heavy RH neutrinos which are singlet under the SM gauge group [44-49] . These singlet 
portal interactions are natural extensions of the SM in the framework of renormalizable quantum field 



^If a Z2 discrete symmetry is the remnant of a broken local symmetry, it can be used to guarantee the stability of dark 
matter as often appears in literature. However the realistic model should contain extra fields and couplings as discussed 
in Ref. [8, 9]. 

^We assume that the local dark symmetry is not confining. The confining nouAbelian hidden sector gauge interaction 
was considered in Ref. [10-15]. 

^ The operators HlhiS are also the SM gauge singlets as H^H, and could be a portal to another singlets from the 
hidden sector. We do not consider this because this operator is dim-5/2, and thus cannot have renormalizable couplings 
with composite operators made of the hidden sector fields charged under symmetry in the sector. Instead, we trade HI Li 
with the lower dim operators Nm's in this paper. 



- 2 - 



theory, and allow rich phenomenology in both dark matter and Higgs sector as we will show in the 
subsequent sections. 

Based on this line of arguments and observations, in this paper we consider an extension of the 
SM where a local U{l)x dark symmetry is introduced to guarantee the stability of dark matter. The 
minimal particle contents and renormalizable interactions are completely fixed once portal interactions 
via Higgs and right-handed (RH) neutrinos are allowed. These extensions allow a possibility to ac- 
commodate neutrino masses and mixings, leptogenesis for BAU, (a)symmetric dark matter and dark 
radiation. In addition to these rich physics, Higgs inflation scenario can be also realized if large non- 
minimal couplings of scalar fields to gravity are introduced, and high enough reheating temperature 
after inflation sets a proper initial condition for the subsequent leptogenesis. 

Before we proceed to the main discussions, let us make two comments. If we considered the 
spontaneously broken U{l)x case by introducing a new {7(l)x-charged scalar cpx with (cpxix)) / 0, 
there would appear a new neutral scalar hx from the radial component of (f)x- Then, this new neutral 
scalar will mix with the SM Higgs h, resulting into 2 neutral Higgs-like scalar bosons. Since hx is the 
SM singlet scalar, the scalar boson sector will be similar to the case of Ref. [42]. There, it was argued 
that the Higgs signal strength is always smaller than 1 independent of the decay channels of the Higgs 
boson, due to the mixing between the SM Higgs and the singlet scalar, and also possible decays of 
scalar bosons into a pair of CDM's. This case would be strongly disfavored if the current situation of 
enhanced — )• 77 remains there in the future analysis. Another issue in the spontaneously broken dark 
symmetry case is the stability or longevity of the dark matter candidate. Nonrenormalizable operators 
suppressed by some powers of (at least) Mp and even renormalizable operators for the scalar dark 
matter case would make the CDM decay in general, as long as electric charge, energy-momentum and 
the total angular momentum are conserved [8, 9]. One has to make a judicious choice of dark charge 
assignments in order to avoid these problems. We postpone the detailed study of the spontaneously 
broken dark symmetry case to the future [9], although we describe the qualitative features in Table 2 
of Sec. 8. 

This paper is organized as follows. In Sec. 2, we define the model Lagrangian assuming the 
local gauge symmetry SU{3)c x SU{2)l x U{\)y x U{1)x as the underlying gauge symmetry, where 
U (l)js: is the unbroken dark symmetry which guarantees stability of the dark matter. The right-handed 
neutrino singlet fields N^iS are also included for the seesaw mechanism and the leptogenesis. In Sec. 3, 
we consider various constraints on our model from large and small scale structure formation, vacuum 
stability and no Landau pole unto Planck scale, direct detection cross section and indirect signatures 
after we identify the dark matter component in our model lagrangian. In Sec. 4, we calculate the amount 
of dark radiation within our model, which originates from massless dark photon. The leptogenesis 
from RH neutrino decays is discussed in Sec. 5. The possibility of Higgs inflation assisted with scalar 
dark matter is discussed in Sec. 6. Collider phenomenology of Higgs boson and scalar dark matter 
is presented in Sec. 7. Some variations of our model is described in Sec. 8, with special emphasis 
on the nature of CDM and singlet portals, the number of Higgs-like neutral scalar bosons, extra dark 
radiation, and the Higgs signal strengths. We discuss a few miscellaneous issues in Sec. 9, including the 
comparison of our model with other models in the literature and effects of nonrenormalizable operators. 
Finally we summarize the results in Sec. 10. Explicit expressions for thermally averaged cross sections 
for the processes relevant to our discussions are presented in Appendix. 

2 The Model 

As explained in Introduction, we assume that dark matter lives in a hidden sector, and it is stable 
due to unbroken local U{l)x dark gauge symmetry. All the SM fields are taken to be U{l)x singlets. 
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Assuming that the RH neutrinos are portals to the hidden sector, we need both a scalar (X) and a 
Dirac fermion (tp) with the same nonzero dark charge (see Table 1). Then the composite operator tpX"^ 
becomes a gauge singlet and thus can couple to the RH neutrinos Njn^s ^. 

With these assumptions, we can write the most general renormalizable Lagrangian as follows: 

^ = ^SM + -^^kin-mix + -^^H-portal + + + -^RHN-portal (2-1) 

where £sM is the standard model Lagrangian and 

/:kin-mix = -^sineX^,^^" (2.2) 

-^H-portal = --^^NxX'^XH'^H (2.3) 

Cx = I (^M + WxqxXf^j - ^X^.X^" - m\X^X - -^\x (^^^)' (2-4) 
= i^jj'' (^d^, + igxqxXf^ ip - m^tpil; (2.5) 



-/:RHN-portai = -MiNgNRi + [Yi^NR^LjH^ + X'NRitl;X^ + H.c.J (2.6) 

gx, Qx, X^ and X^i, are the gauge coupling, U {l)x charge, the gauge field and the field strength tensor 
of the dark U{l)x, respectively. B^i, is the gauge field strength of the SM C/(l)y. We assume 

> 0, Ax > 0, Xhx > 0, (2.7) 

so that the local U{l)x remains unbroken and the scalar potential is bounded from below at tree level 

5 

Either X or is absolutely stable due to the unbroken local U{l)x gauge symmetry, and will be 
responsible for the present relic density of nonbaryonic CDM. In our model, there is a massless dark 
photon which couples to the SM U{1)y gauge field by kinetic mixing. One can diagonalize the kinetic 
terms by taking a linear transformation defined as [28] 




1 / cos e 
— tane 1 




(2. 



In this basis, the SM U{1)y charge is redefined as qy = qy/ cose, and hidden photon does not couple 
the SM fields. However, dark sector fields now couple to the SM photon and Z-boson. In the small 
mixing limit, the couplings are approximated to 



>Cds-sm = ipij'' [d^ - igxqxte [cw^^i - swZf,)] V' + - igxqxte {cwA^ - swZ^)] X\ (2.9) 

where = tan e, cw = cos 9w and sw = sin 9w with Ow being the Weinberg angle. Hence, dark sector 
fields charged under U (l)x can be regarded as mini-charged particles under electromagnetism after 
the kinetic mixing term is removed by a field redefinition, Eq. (2.8). 

Meanwhile, we can assign lepton number and U (l)x charge to RH neutrinos and dark fields as 
shown in Table 1. Then, the global lepton number is explicitly broken by Majorana mass terms for the 
RH neutrinos. If Yjj and Aj carry CP-violating phases, the decay of RH neutrinos can develop lepton 



''if we did not assume that the RH neutrinos are portals to the dark sector, we did not have to introduce both ^ and 
X in the dark sector. This case is discussed in brief in Sec. 8. 



^Quantum corrections to the scalar potential will be discussed in Sec. 3.2. 
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field 


N 




X 




1 


1 










1 


1 



Table 1. Lepton number and U{l)x charge assignment 

number asymmetry in both of visible and dark sectors. Since U{l)x is unbroken, the asymmetry in 
the dark sector has a relation, 

Ya^ + Yax = (2.10) 

where Yai = {ni — n^) /s is the asymmetry between i and i with rii and s being the number density of 
i and entropy density. 

There are various physics issues involved in our model as listed below: 

• Small and large scale structure 

• Vacuum stability of Higgs potential 

• CDM relic density and direct /indirect DM searches 

• Dark radiation 

• Leptogenesis 

• Higgs inflation in case of a large non-minimal gravitational couplings 

In other words, the model will be highly constrained, but astonishingly it turns out that our model 
can also explain various issues related to those physics in its highly constrained narrow parameter 
space without conflicting with any phenomenological, astrophysical and cosmological observations. It 
is highly nontrivial that our model can accommodate all these constraints in a certain parameter 
region, reminding us that our model was based on local gauge principle for the dark matter stability, 
and assumption of singlet portals to the dark sector, by introducing only 3 new fields, X,ip and X^. 

3 Constraints 

Including the portal interactions, the presence of an unbroken local U{l)x in dark sector with kinetic 
mixing with the SM sector is subject to various phenomenological and cosmological constraints. In 
this section, we will take a look each of constraint or physics one by one. 

3.1 Structure formation 

The presence of the dark matter self-interaction caused by nonzero charge of U {l)x could affect sig- 
nificantly the kinematics, shape and density profile of dark matter halo, so it is constrained by, for 
example, the galactic dynamics [50], ellipticity of dark matter halos [51] and Bullet Cluster [52] (see 
also [19, 25, 53, 54]). For a velocity-dependent self-interaction, the transfer cross section of the dark 
matter self-interaction, defined as ctt = / dfl (1 — cos9) is upper-bounded as [25] 

< 35 cmVg (3.1) 

ti=10km/s 



a, 



obs 
T 
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Interestingly, it is shown that, if a^j^'^ is close to the bound, it can solve the core/cusp problem [21] and 
"too big to fail" problem [24] of the standard collisionless CDM scenario [25]. 

In our model, for both -0 and X the self-interaction cross section with a massless dark photon is 
given by [55] 



(47rpxa|)i/2 



(3.2) 



where v and px are the velocity and density of the dark matter at the region of interest ^. We take 
V = lOkm/sec and px = 3GeV /cm^. Then, compared to Eq. (3.1), dark interaction is constrained as 

where we approximated the log factor to 41. Either X or ■0, which is lighter than the other, poses a 
stronger constraint on ax- Note that tp couples only to dark photon at low energy, and the thermally- 
averaged annihilation cross section of is found to be 

The abundance of ip at freeze-out is 



" T,,^ ' [ 9*s J \Ti,^J {av)tnnMp " {avf ^ " ' ' 
where 

(f^^')ann ~ 6 X lO^^^cm^sec. (3.6) 

is the thermally-averaged annihilation cross section which gives the right amount of present dark matter 
relic density ^ corresponding to 

~ 2 X 10-10 GeV (3.7) 

S / obs 

In the far right-hand side of Eq. (3.5), we used the fact that , even if {av)tnn varies by several orders 
of magnitude, mdm/^f is changed only by a factor of 0(1). The constraint, Eq. (3.3), implies that 

{crv)tnn [ (^j for < mx 

Hence, if were stable, it would be over-abundant at present. 
In order to avoid the over-closing by ■0) we assume 

> mx (3.9) 

so that ijj can decay through the virtual RH neutrinos. The decay rate of ip is given by 

— T^^uX + '^ip^uXh (3.10) 



ann 



obs 



''There are other t-channel scatterings of X-X'' (the contact interaction- Ax, Higgs and Z-boson mediations), but they 
don't have large enhancement caused by small velocity. 

^ Since a dark matter charged under an unbroken symmetry can annihilate only with its anti-particle which constitutes 
the half of the whole CDM relic density, ((7t;)ann is larger than the one for a charge-neutral dark matter by a factor 2. 
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where 





) 



2 




167r Ml 



(3.11) 




487r2 I V- 



(3.12) 



with fhu = Y^vjj/Mi and vh = 174 GeV being respectively a contribution to the neutrino mass matrix 
and the vev of Higgs field. The present CDM relic density poses the strongest constraint on F^, as we 
win see in a moment. 

Eq. (2.10) imphes that, even if the asymmetry between and ip may arise in the decay of RH 
neutrinos, once ip decays, the dark matter composed of X and becomes totahy symmetric irrespec- 
tive of its origin. If ip decays before the thermal component of X freezes out, the X's coming from 
the decay of ip thermalize, which makes the number density nx return to that of thermal equilibrium. 
The present relic density in this case is determined by the thermal relic, hence the annihilation cross 
section should be the one in Eq. (3.6) ("symmetric thermal" case). On the other hand, if '0 decays after 
the thermal freeze-out of X, the annihilation cross section should be larger than the one for thermal 
relics (i.e., Eq. (3.6)) so that the non-thermal freeze-out to provide a right amount of relic density 
("symmetric non-thermal" case). In this case, the required background temperature when -0 decays is 
determined by the annihilation cross section of X. 

In our model, the pair annihilation of X-X^ can be controlled by the Higgs portal interaction Xhx 
which leads s-wave annihilations. It freezes out at a temperature Tf ~ mx/20. However dark matter 
can still be in kinetic equilibrium with thermal background at a lower temperature due to scatterings 
to SM particles. The scattering is mediated by photon and Higgs thanks to the kinetic mixing and 
Higgs portal interaction. The transfer cross section of photon- mediation is such that ax oc e'^/T'^. 
Although the associated scattering could be quite efficient at a low temperature, e ^ 1 make it less 
efficient. We found that, for e ~ 0(10"^) which will be of our interest as described in section 3.3, the 
momentum transfer rate via photon is too small to keep kinetic equilibrium after freeze-out. In case of 
Higgs mediation, the kinetic equilibrium can be maintained by the scattering mainly to charm quark 
to a temperature of the charm quark mass scale [56]. At a lower temperature, the scattering rate is 
too small. 

Hence, for Xhx ^ 1 and electroweak scale mx, the kinetic decoupling takes place before QCD 
phase transition ^. If ip decays to X abundantly, X and would be able to re-annihilate even after 
freeze-out of the thermal annihilation until their number densities is reduced enough to stop the re- 
annihilation. The abundance of X and X^ at the moment should be responsible for the present relic 
density of dark matter. Hence, when ^l; decays at a temperature T^, the annihilation of X should be 
frozen with a rate. 



being the present number density to entropy density ratio of X that matches observation, (o'w)^^^^ is 
the annihilation cross section of X when ip decays. Equating the annihilation rate to the expansion 

* As long as X is decoupled before QCD-phase transition, the effect of the dark photon to the SM radiation at the time 
of BBN is negligible even though the dark photon is decoupled from X at temperature T ~ 16MeV [^^^ — ) ( ano&iv )^'^^ 




(3.13) 



with 




(3.14) 




in our scenario [55]. 
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rate when ^/^ decay, we find the decay temperature of ijj to give a right amount of non-thermal rehc 
density, 



Th = 



-1/4 



V^r^Mp ~ 9.2 GeV 



th 
ann 
X 



(3.15) 



where we used fif*(rd) = g^siTd) = 100 and Mp = 2.4 x 10^*^ GeV in the right-hand side of above 
equation. This implies that 

'0.1eV\ / Ml \ / i lev \ m 



Af ~ 58.5 




X 



1 + 



1 



487r^ V^J? / 



300 GeV 



(3.16) 



annd ~ ('^^)anni equal to or larger than Tf and corresponding Ai are fine. Note 



Note that, if {av)^ 

also that > T^d unless {cv)^^^ ^ is larger than(iTz;)*nn by at least two orders of magnitude. However, 
as described in section 3.3, only (o'w)^n/(fTt;)!;^^ ^ 5 is allowed, so we can take {o'v)^^^^ = {(tv)-^^^. 
Fig. 1 shows contours for a right amount of relic density as a function of Ai and m^. In the figure, 
solid and dashed blue lines are respectively for (o'w)ann/(o"w)ann = ^ 1 with = 0.1 eV, Mi = 
1.63 X 10^'^ GeV and nix = 300 GeV. We also depicted decay temperature of tp there. 

Shortly speaking, the existence of the massless dark photon constrains our model parameters to 
satisfy Eqs. (3.3) and (3.16). They are from small/large scale structure formation and present dark 
matter relic density, respectively. 

3.2 Vacuum stability 

In the standard model, Higgs potential becomes unstable at an intermediate scale because of top 
loop contributions to the Higgs quartic couplings, though it depends on some of the standard model 
parameters, for example top pole mass and strong interaction [57]. Such instability can be cured if 
Higgs field couples to other scalar field(s) [58-60] . Depending on the existence of mixing between Higgs 
and additional scalar (s), tree-level and/or loop effects should be able to remove the instability. In our 
model, X does not develop non-zero VEV, and the SM Higgs is not mixed with X. In this case, the 
loop-effect should be large enough to remove the vacuum instability of the SM Higgs potential. Note 
that the Dirac neutrino mass terms also contribute to the RG-running of the Higgs quartic coupling. 
However it is a negative contribution reflecting the fermionic nature of the right-handed neutrinos [61]. 
Hence, in order not to make worse the vacuum instability up to Planck scale, we take 



y:' < 0.1 



(3.17) 



and ignore its contribution to the RG equation of Higgs quartic coupling. Then, the relevant one-loop 
RG equations are 



dln/i 



where i = H, HX, X and 



/3. 
/3a 



1 

^" ~ 16^ 

_ Affx 
~ 167r2 

1 



2A\l + 12AhA? - 6Af - 3Ah {2,gl + gl) + ^ {2gl + (gl+glf) + ^Xjjx 



2 (6Ah + 3Ax + Xhx) - i:Xh {^gl + gl) - 6A: 



167r2 



2'^HX + 18Ax 



(3.18) 

(3.19) 
(3.20) 
(3.21) 
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mx=300GeV, mv=0.1eV, Mi = 1.63xlO'^GeV 




m^[TeV] 

Figure 1. Parameter space for a right amount of dark matter relic density. Contours correspond to the present 
dark matter rehc density for (o'w)^j„/(crw)^^j,j = 1,5 (sold blue lines from right to left) with nix — 300 GeV, 
nil, = 0.1 eV and Mi = 1.63 x 10^" GeV. The gray region is excluded by XENONIOO direct dark matter search 
as described in section 3.3. 

in addition to the ones for the other SM couplings. We solved 2-loop RGEs for SM couplings and 
1-loop RGEs for non-SM couplings numerically, and found that the vacuum stability of Higgs potential 
and perturbativity of the couplings requires 

0.2 < Xhx < 0.6, Xx < 0.2 (3.22) 

3.3 Direct detection 

In our model, dark matter couples to the SM particles via neutral SM gauge bosons (see Eq. (2.9)) 
and Higgs portal, hence both type of interactions provide channels for dark matter direct searches. 
In the case of gauge boson exchange, the spin-independent (SI) dark matter-nucleon scattering cross 
section via photon exchange provides a strong constraints on the kinetic mixing. As can be seen from 
Eq. (2.9), our dark matter has a mini-electric charge, 

Eg = — — gxcpi/ tan e (3.23) 
e 
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For a scattering to a target atom with atomic number Z, the differential cross section of the Rutherford 
scattering of our dark matter is given by 

ail A:m\v^^suC{0 12) 

where = T^xf^A/ {mx + "ia) with ruA being the mass of the atom is the reduced mass, is the 
dark matter velocity at the center mass frame, and J-A{<irA) is the form factor of the target atom with 
q and being respectively the momentum transfer and effective nuclear radius. The scattering angle, 
0, is related to the nuclear recoil energy of the atom, as 

= — t^^(l-cos0) (3.25) 
ruA 

where v is the lab velocity. So, Eq. (3.24) is expressed as 



da A _ 2T:elal^Z^ ^ 



TX{Er) (3.26) 



For the SI dark matter-nucleus scattering, the differential cross section with respect to the nucleus 
recoil energy is parameterized as 

(3.27) 

where Hp = mxrup/ {mx + rup) is the reduced mass of dark matter-proton system, (cp^)g^p is the 
dark matter-proton scattering cross section constrained by experiments. Hence, the kinetic mixing is 
constrained as 



te < 



1 



1/2 



where EJ is the threshold recoil energy of a target atom at a given experiment. 
In the case of the Higgs portal interaction, the scattering cross section is 



<^U,H = ^^rfkn (3-29) 



where 



fj\f,H = z^HX ^-^fg,H (3.30) 

8 mxmjj 



with 

fq,H 



q=u,d,s q=t,b,c 



(3.31) 



and and /q being hadronic matrix elements with a scalar. Based on the study on lattice [62], we 
take fq = 0.326 here. Currently, the strongest bound on comes from XENONIOO direct search 
experiment [63] which has EJ = 6.6 keV. Fig. 2 shows how the kinetic mixing (left panel) and 
Higgs portal coupling (right panel) are limited by the experiment (gray region). Also, depicted are 
thermally averaged annihilation cross sections (colored lines) and bounds from vacuum stability and 
perturbativity (dark gray regions). In the left panel of the Fig. 2, we notice that, if small scale anomalies 
of structure formation are to be explained by the dark matter self-interaction, XENONIOO direct search 
experiment constrains strongly the kinetic mixing as 

e < 10"^ - 10-^ for 6 GeV < mx < I TeV . (3.32) 
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10' 10- I0-' 10' 10- 10' 



Figure 2. Left: XENON100(2012) bound on the kinetic mixing parameter e as a function of mx for ax given 
by the bound value of Eq. (3.3). Right: XENONIOO (2012) bound and contours of {crv)^^^J {av)f^^ in {Xhx, mx) 
plane. ((Tw)ann — 2 X 10~'^^cm^ is the thermally averaged annihilation cross section giving the correct amount 
of dark matter relic density from thermal freeze-out. The gray region is excluded by the recent result from 
XENONIOO [63]. The lower and upper dark gray region is excluded by the vacuum stability of Higgs potential 
and perturbativity of couplings, respectively. The colored lines correspond to (o'w)^in/ (o'w)ann — 1' 2, 5, 10 from 
bottom to top. 



From the right panel of Fig. 2, we also notice that direct search experiments already excluded mx ^ 
80 GeV except for the narrow resonance band around mx — 'mh/2. In addition, for mx = 0(10^^^) GeV, 
(fjf can be larger than (the one for the right amount of thermal relic) only by an order 

of magnitude at most. For mh = 125 GeV, if top pole mass and strong coupling are respectively 
mt = 173.2 GeV and = 0.1184, vacuum stability and perturbativity allows mx only in the range 

200 GeV <mx< 600 GeV (3.33) 

and annihilation cross section satisfying 



\ /ann 



< 5 (3.34) 



This implies that the thermal relic can be reduced to abut 20 % of the present relic density at most, 
and asymmetrically produced non-thermal dark matter can saturate the present relic density. Note 
that the recent report on ~ 130 GeV line spectrum in Fermi-LAT 7-ray data is not achievable in 
our model since the branching fraction dark matter annihilation to photon(s) is of O(10~^ — 10~^). 

3.4 Indirect Signatures 

The dark interaction and kinetic mixing in our model should be highly suppressed as described in 
previous sections. In addition, since ax ^ 10"^ for mx ^ 0(1) TeV (see Eq. (3.3)), Sommerfeld 
enhancement factor, which is given by 

S = ^"^/^ (3 35) 

1 — e~'^'^xh ' 



- 11 - 



is 0(1). Hence it is difficult to expect detectable indirect signatures from the annihilation channels via 
dark interaction or kinetic mixing. 

The possible indirect detection signatures comes from Higgs portal interactions, 

XX^ ^ H* ^ ff, VV, or XX^ HH, 

where / and V are the SM fermions and the weak gauge bosons, respectively. These processes can 
produce a sizable continuum spectrum of photons, since the annihilation cross section can be larger than 
the value for thermal dark matter. However, the recent data from Fermi LAT 7-ray search provides 
upper-bounds on various annihilation channels [64]. In our model, W~^W~ channel is dominant. Taking 
into account the fact that an annihilation is possible only for X-X'^ pairs, the annihilation cross section 
is expected to be constrained at least as [64] 

Mxxt^VK+H/- < 2 X 7.4 X 10-26cm3/sec (3.36) 

for NFW dark matter profile. Hence the total annihilation cross section is upper-bounded as 

{'^v)^nn ^ Br(XXt ^ W+W-y^ X 2 X 7.4 X lO^^^cm^sec (3.37) 

In the allowed region of parameter space, that is, for mx = 0(10^^^) GeV, we find Br(XX^ — ?■ 
VF+iy^) ~ 0.5, and the allowed ratio of the annihilation cross section to the value for thermal relic is 
bounded as 

1 < < 5. (3.38) 

This constraint is similar to the one coming from the perturbative bound shown in Fig. 2. 



4 Dark Radiation 

Dark photon can contribute to the radiation density of the present universe. Its contribution is pa- 
rameterized in terms of the extra relativistic neutrino species as 

AiVeff = ^ = ^ f^V i ^#4V'' (4-1) 



Tufl) \T^i, dec J \g*s{T^Acc 
where and pj/ are respectively the present energy densities of the dark photon and a neutrino species, 
Qi, Tifl and Tj dec are respectively the degrees of freedom, the temperature at present and decoupling 
of the species, i, and g^s is the total SM degrees of freedom associated with entropy. Because of the 
energy injection to photons at the epoch of electron-positron pair annihilation which took place after 
neutrino decoupling, the photon is slightly hotter than neutrinos at present, resulting in the ratio of 
temperatures, T^^q/T^a = (4/11)^''^. In addition, dark matter is decoupled from the SM thermal bath 
at a temperature T ~ 1 GeV before QCD-phase transition while still in contact with dark photon. 
Hence dark matter and dark photon are decoupled from the SM thermal bath simultaneously. When 
it is decoupled, the temperature of dark photon is the same as that of photon. Therefore, we find 

AiVeff = f ^^f^i°\y^' ^ 0.06 (4.2) 

where we used g*s{T^fl) = 3.9 and 5*5(T^,doc) = 100. The best fit value of observations is [26] 

Ar°|s = 3.84 ± 0.40 at 68% CL (4.3) 

with SM expectation = 3.04. Therefore, in our model the contribution of dark photon to the 

radiation density at present is consistent with observation within about 2-cr error, slightly improving 
the SM prediction in the right direction. 
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5 Leptogenesis 



Our model allows production of lepton number asymmetries in both of visible and dark sectors via 
decays of heavy RH neutrinos. If the mass of dark matter X is much larger than proton mass and 
asymmetric generation of dark matter is responsible for the present relic density, the asymmetry of ip 
should be much smaller than that of lepton £j. However, the contribution to X and from the decay 
of thermal symmetric component of ip-ip is dominant as described in section 3.1. The present relic 
density is then determined by thermal or non-thermal freeze-out of the annihilation cross section of 
X, depending on the temperature when tp decays. Considering asymmetric generation of dark matter 
in this circumstance is pointless. However we still have to check if a right amount of lepton number 
asymmetry in the visible sector can be achieved. 

The lepton number and U{1) charges are assigned to relevant fields as shown in Table 1. Then, the 
global lepton number is explicitly broken by Majorana mass terms for the RH neutrinos. The lightest 
RH Majorana neutrino A^i can decay into both the SM fields and the DM fields: 

With nonzero complex phases in Yy and Aj the decay can generate the AL, and /S.X as ^ 

^A. ^ ^ = e,mY^\Q) (5.1) 

where is the number density of a charge asymmetry associated with the field i, s is the entropy 
density, and r]iS are asymmetry and wash-out effect of field i from the decay of A^i, respectively, and 
y^^^(O) = 135C(3)/47r'^5r* with g^{T > Mi) ~ 100 being the number of relativistic de grees of freedom 
at a temperature well above the mass scale of the lightest RH neutrino (Mi). For hierarchical mass 
spectrum, Mi <C M2,3, the asymmetries are given by [49] 



Ml I™ 


(3y;yJ + a*a) m 


-'y^yJ 


11 


Svr 


2YyYj + AA* 


11 



Ml Im [{Y;yJ + A* A) M-iAA*] 



Svr 



^'^ 8^ \2YyYj + AA* 

L J 11 

where M = diag (Mi, M2, M3), and upper-bounded as [49, 65] 



(5.2) 
(5.3) 



^ 3Mim^^^ ^ ( 1 for Br^ > Br^ 

- WTTvj^ ^ [ ^Ap'fi/A^Ms for Brz, < Br^ ^ ' ' 

with vh = 174 GeV being the VEV of the SM Higgs and m™^^ is the mass of the heaviest left-handed 
neutrino. 

The visible sector lepton number asymmetry, AL, would end up the visible sector baryon number 
asymmetry via anomalous electroweak process [66, 67]. For simplicity, if we assume the visible sector 
lepton number asymmetry is dominated by a flavor, the late-time baryon number asymmetry is related 
to the lepton number asymmetry as [2] 

Yab = ^Ial (5.5) 
Then, the present observations of baryon number asymmetry can be matched if 

Yal ^ 2.6 X 10"^° (5.6) 



For simplicity, we do not consider the case where the initial abundance of A^i is negligible or zero. 
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The eventual outcome of leptogenesis via the decay of heavy RHN can be obtained by solving 
Boltzmann equations which involve effects of wash-out and transfer of the asymmetries between visible 
and dark sectors. However, if the narrow-width approximation is valid, we can get much simpler 
picture. The narrow-width approximation is valid if 

\r « 1 (5-7) 



MiHi 



and 2 — )■ 2 scattering between visible and dark sectors via heavy neutrino is ineffective, hence asym- 
metries in both sector evolves independently. In this circumstance, the washout effect on asymmetry 
is mainly from the inverse decay. If the washout effect is weak, i.e., 

< 1 (5.8) 

-Hi 

the final asymmetry is directly related to the asymmetry from the decay of RHN. Otherwise, there can 
be large reduction of the asymmetry. 
The decay rate of RHN is 

Ti = (n'l + A?) Ml (5.9) 



and the branching fractions to the SM and dark sectors are 



^'^i Br , - ^ 



Hence 



£]_ _ Mp ( my/vjj for L 
'iri~\M^\ Af/Mi for V 



Br.^ = 9^ X <! ' z : (5-11) 



For simplicity, we use narrow-width approximation from now on. This implies 



Note that 



= ^ - 3 X 10- ( ) ( ) (5.13) 



rfivM-i „ a ( "fnv \ ( M- 



vi 



O.leV; VlO^GeV 



So, Yy\ can not saturate the bound of Eq. (5.12) for M\ <C 10^^ GeV which we assumed in order not 
to worsen the vacuum instability of the SM Higgs potential. Hence, for Ai saturating the bound, we 
always have Br/, <C Br^. Combined with the constraint Eq. (3.16), the narrow- width approximation 
can be achieved if 



> 94.3 TeV 



my ) VlO^GeV ) \{av) 



0.1 eVW Ml \'/V(^^)ann 



1/3 

) (5.14) 
300 GeV ' ^ ' 



Depending on the sizes of Yy\ and Ai, there are various regimes of wash-out as analyzed in Ref. [49]. 
The purpose of this paper is not at the full analysis of leptogenesis, so here we simply show a working 
example in the following paragraph. 

If ihy ~ 0.1 eV and Ai > Yv\^ both of visible and hidden sectors are in the strong washout regime. 
The wash-out effects are given by [2, 49] 

TiBrz, TiBr^ 
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with the ratio between asymmetries, 



(5.16) 



Since can be smaller or larger than Ial even though mx is much larger than proton mass, we 
can assume 

^ = 1 (5.17) 
From Eqs. (5.1), (5.4) and (5.15), the maximally expected late-time lepton number asymmetry is 

Hence the present baryon number asymmetry corresponding to Eq. (5.6) can be obtained if 



where we used Eq. (5.17) in the left-hand side of above equation. Fig. 3 shows a parameter space 
limited in our analysis. In the figure, darker gray region is excluded by XENONIOO direct dark matter 
search. Narrow-width approximation is valid in the white region well below the border line of the 
lighter gray region (e.g., below the dashed gray line). Ai < Yi^i below the green line. Although a 
wider parameter space may be allowed, our analysis of leptogenesis in this section is limited only in the 
white region bounded by the border of the lighter gray and green lines. In the region, right amounts 
of baryon number asymmetry and dark matter relic density can be obtained as long as Eq. (5.19) is 
satisfied. 

So far, we have considered the lepton number asymmetry in the visible sector that comes from the 
decay of RH-neutrinos only. However there is an additional contribution from the late-time decay of 
^ which also carries lepton number. Since the decays of ■0 and ip involve a virtual internal line of a 
Majorana RH-neutrino which decays eventually to a SM lepton and Higgs pair, both of decays produce 
equal amount of the same-sign lepton number asymmetry in the visible sector. In addition, there is no 
dilution of the produced visible sector asymmetry due to inverse decay or transfer to the dark sector, 
since such processes are kinematically forbidden. Hence, the contribution from those decays is 

A{Yal) = 2eLY^(Tfj (5.20) 

where is the freeze-out temperature of the pair annihilation of V'-V' and we used Y^{T^) = Yjp{T^) 
in the right-hand side of the above equation. The freeze-out abundance of is given by [68] 



= 3.79(V8^)-^yVg..xt ^ 0.054?? (5-21) 



where xf^ = m^/T^ and we used g^, = g^:S = 100 and Eq. (3.4) at the far right-hand side of the above 
equation. Combining with Eqs. (5.4) and (5.6), we find 



A(1al) ^ 2 X ^nl^m^ M^mr^ ^ I 1 ^ for Br^ » Br^, ^2 

Yal a\ Mp v]j \ y/XiMi/XfA'h for Brz. < Br ^ " 



As an example, we may take ei ~ 10"''. Then, for ax = 10~^ and = 10^ TeV, we find xf^ ~ 2.2 
resulting in A(1al)/^al — 0.3. Therefore, depending on ax and m^, the decay of ip and -ip can be the 
origin of the baryon number asymmetry in the present universe even though the asymmetry between 
-0 and tp is absent. 
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mx=300GeV, mv=0.1eV, Mi = 1.63xlO'^GeV 




10° IQi 10^ 10^ 

m^[TeV] 

Figure 3. Parameter space for right amounts of baryon number asymmetry and dark matter relic density at 
present. We used mx = 300 GeV, = 0.1 eV and ^/Y^M^Jy^^ = 1 corresponding to Mi = 1.63x101" GeV. 
Dark gray region is excluded by XENONIOO dark matter direct search experiment. In the dark gray region, 
narrow-width approximation is not valid. The boarder of the light dark gray region and the gray dashed line 

correspond to Ai/ \J IQtt^ Mi / Alp — 1, 1/3, respectively. Below the green line, Y^i > Ai for which our analysis 
is not valid. The blue lines correspond to (o'w);^n/(^^)ann = from right to left. 

6 Higgs Inflation 

In order for the leptogenesis described in the previous section to work, the temperature of the early 
universe should be high enough so that the lightest RHN can be in thermal equilibrium before it is 
decoupled. This condition can be achieved if the reheating temperature of the primordial inflation is 
high enough. An intriguing possibility is so-called Higgs inflation [69, 70] which uses the SM Higgs as 
the inflaton equipped with a large non-minimal gravitational coupling. As a variant, Higgs-scalar singlet 
system has been also considered in the literature [71] (see also [72, 73]). Modulo the subtle issues of the 
unitarity problem [74], our model indeed allows inflation along Higgs direction since Higgs potential is 
stabilized by the help of a coupling to the singlet scalar X. The model parameters relevant to inflation 
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are Xhx, and the Higgs quartic coupling in addition to the large non-minimal couplings (say ^j). 
As free parameters, we can adjust ^jS for given set of quartic couplings while satisfying requirements 
on the inflationary observables under the assumption of the positivity of quartic couplings (see [58] for 
example). Hence the physics involved in inflation does not pose any new constraint other than ones 
described in previous sections if inflation takes place along Higgs direction, and the Higgs inflation 
along with a singlet scalar can be realized. 

It turned out that the reheating temperature after Higgs inflation is around 0(10^'^"^^) GeV [75]. 
It is high enough to populate the lightest RHN in thermal bath. Therefore, Higgs inflation sets the 
initial condition for the leptogenesis. 

7 Higgs and DM phenomenology at colliders 

The Higgs boson in our model could decay into a pair of scalar DM's through Xhx term if kinematically 
allowed. However, as shown in Fig. 2, dark matter direct search allows only mx ~ m]^j2 with Xhx ^ 
10~^ even though SM Higgs may not suffer from vacuum instability problem. If it is allowed, the decay 
rate of Higgs to dark matter is 



and we found that the branching fraction of Higgs to dark matter can be consistent with recent LHC 
data only if Xhx ^ 0.1 or mx is very close to mij2. On the other hand, if vacuum stability is imposed, 
such a small Xhx is excluded and only mx = 0(10^"^) GeV is allowed. In this case, the production 
and decay rate of Higgs boson in our model are exactly the same as those of SM Higgs boson, since 
H — )■ XX^ is kinematically forbidden. Therefore it is difficult to discriminate our model from SM. In 
other words, if collider experiments shows any non-SM signature, our model is excluded. 



It may be possible to search for a pair of dark matter production at the LHC or the ILC through 
e+e" ^ Zh* Z{XX^), or WW fusion through 



with extra emissions of gluon or 7 from the initial or the final quark jets. The detailed study of this 
channel will be beyond the scope of this paper, and will be addressed elsewhere. 

8 Variations of the model 

Instead of our model analyzed in this paper, one can consider a simpler dark sector which contain either 
X or ^ only in addition to X^. In these cases, renormalizable RH- neutrino portal interactions are not 
possible and leptogenesis from seesaw sector has nothing to do with the dark matter. If ip is absent and 
X, X^ are dark matters, the only change relative to our present model is that the current dark matter 
relic density should come from the thermal freeze-out of X-X^ annihilation via Xhx interaction. Hence 
the annihilation is fixed to be Eq. (3.6) as usual. 

If X is absent and ip is the dark matter, one has to introduce a real SM-singlet scalar (say 5) 
connecting the dark sector to the SM sector as the model discussed in Ref. [42], so that the thermal 
freeze-out of ip-ijj annihilation via the newly introduced interactions of S provides a right amount of 
dark matter at present. Otherwise ip and would be overproduced due to the smallness of ax- The 
physics of this model is nearly same as that discussed in Refs. [42, 60] modulo the effects of the dark 
interaction on structure formation and direct DM searches, as well as dark radiation from massless 




(7.1) 



qq qqh* qqXX^ 
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hidden photon. The spin-independent cross section of the ip (or ^)-to-nucleon scattering via photon 
exchange is the same as the one in the case of X-X'^ dark matter, so the constraint on the kinetic 
mixing (e) shown in the left panel of Fig. 2 is equaUy applicable to this case. Higgs inflation is still 
possible in these variations as discussed in Ref. [58], since there are extra scalar fields X or S" in either 
case. 

Finally one could consider the case the U{l)x dark symmetry is spontaneously broken by nonzero 
{<i>x) 7^ 0. Then in this case there is a singlet scalar from (px after U{l)x breaking, which will mix 
with the SM Higgs boson. Therefore there are two Higgs-like neutral scalar bosons after all, and both 
of them have signal strengths universally suppressed from the SM value "1" . If tpx is the CDM, one 
needs a singlet scalar 5 as a messenger, and this will mix with the SM Higgs boson (and the remnant 
from 4)x)- 

In Table 2, we summarize the dark field contents, messengers, the particle identity of the dark 
matter (DM), the amount of dark radiation (DR) and the signal strengths of Higgslike neutral scalar 
bosons (including the number of them) in various scenarios. In all cases, there are additional scalar 
bosons (either X or (px or both) which make Higgs inflation still viable for mn = 125 GeV. And the 
Higgs signal strength is smaller than "one" except for the scalar is the CDM with unbroken U{l)x dark 
symmetry. Especially /^i=i,2,(3) < 1 for fermion CDM, whether U{l)x is broken or not. Our conclusions 
on the Higgs signal strength is based on the assumption that there is only one Higgs doublet in the 
model. If we include additional Higgs doublets or triplet Higgs, the Higgs portal would have richer 
structure, and the signal strength will change completely and will vary depending on the Higgs decay 
channels. Also it should be possible to have a signal strength for H ^ channel greater than "one" 
without difficulty. 



Dark sector fields 


U{l)x 


Messenger 


DM 


Extra DR 




Xf„X,iJjx 


Unbroken 


H^H,B^,X^''',Nn 


X 


~ 0.06 


1 (i = 1) 


Xfi, X 


Unbroken 


H^H, B^,X^'^ 


X 


~ 0.06 


1 (i = 1) 


X^,ipx 


Unbroken 


H^H,B^,Xf''',S 


ipx 


~ 0.06 


< 1 (i = 1,2) 


X^,X,7px,(px 


Broken 


WH,B^,Xf'\NR 


X or ^px 


~ 


< 1 (i = 1,2) 


Xfj,,x, 4>x 


Broken 


H^H,B^,X^'^ 


X 


~ 


< 1 (i = 1,2) 


X^,ipx 


Broken 


H^H^B^.X^'^S 


ipx 


~ 


< 1 (i = 1,2,3) 



Table 2. Dark fields in the hidden sector, messengers, dark matter (DM), the amount of dark radiation (DR), 
and the signal strcngth(s) of the i scalar boson(s) (/^i) for unbroken or spontaneously broken (by {4>x) i^ 0) U (l)x 
models considered in this work. The number of Higgs-like neutral scalar bosons could be 1,2 or 3, depending on 
the scenarios. 

9 Discussions of some miscellaneous issues 
9.1 Comparison with other models 

Leptogenesis in our model is very similar to one in Ref. [49] in that the RH neutrino decay is the origin 
of BAU and CDM observed today. However there are a few important different aspects of our model 
compared with Ref. [49]: 

• Our lagrangian is based on local gauge symmetry Gx = U{l)x that guarantees that CDM is 
absolutely stable. Assuming all the SM singlet fields are portals to the hidden sector DM, we are 
naturally led to the present model without any other ad hoc assumptions. 
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• If stable, either X ov ^ could be dark matter, but the smallness of the dark interaction that is 
required from observations of large scale structure does not allow the fermion dark matter ip in 
our model. Hence ip should be able to decay. The interaction to SM via RHN allows the decay of 
tp ii X is lighter than ip. Due to this process, dark matter is composed of X and X^ , and becomes 
symmetric eventually irrespective of its origin (symmetric thermal or asymmetric non-thermal). 

• Because of the smallness of dark interaction, the freeze-out abundance of ip before decay is 
quite large and becomes main contribution to dark matter abundance of X-X^ . In other words, 
asymmetric production of dark matter does not play any significant role in our scenario, and the 
eventual relic density of dark matter is determined by thermal or non-thermal freeze-out of X-X^ 
pair annihilation through Higgs portal XxH terms. 

• The decays of ip and ip via RH neutrino portal contribute to the visible sector lepton number 
asymmetry even if there is no asymmetry between ip and ip. The contribution can be the origin 
of the present baryon number asymmetry if the coupling of the dark U{l)x is small {ax ^ 10^^) 
and V' is heavy enough (m^ > 10^ TeV). 

• Higgs inflation can be realized thanks to the existence of the portal interaction which is necessary 
for efficient pair annihilation of dark matter X and X^ and for vacuum stability. The large enough 
reheating temperature after inffation can set a proper initial condition for leptogenesis to work. 

• Since the dark symmetry is unbroken gauge symmetry in our model, there is always dark radiation 
from massless hidden photon. This conclusion can be evaded, if the dark symmetry is an unbroken 
but confining symmetry like color gauge symmetry in QCD [10-13]. In that case, CDM would 
be a composite hadrons made of hidden sector quarks (similar to baryons or mesons in ordinary 
QCD without electroweak interacrtion) and would be absolutely stable. 

• The dark matter self-interaction caused by the massless dark photon can explain the small scale 
structure problem appearing in usual collisionless CDM scenario. 

9.2 Effects of nonrenormalizable operators 

Since our model has no Landau pole or vacuum instability up to Planck scale (Mpi), this model could be 
an ultimate theory up to Planck scale when we ignore the fine tuning problems related with the Higgs 
mass^ or cosmological constant. Still there may be higher dimensional nonrenormalizable operators 
suppressed by some positive powers of 1/Mp originating from quantum gravity effects. Since the dark 
symmetry U{l)x is not broken, dark matter would be absolutely stable even in the presence of these 
higher dimensional operators 

In this section, let us list some dim-5 or dim-6 operators that are suppressed by one or two powers 
of 1/Mp and contain either the RH neutrino or the dark fields and discuss their effects on the results 
obtained in the previous sections. We still impose the local gauge invariance of the nonrenormalizable 
operators under local gauge symmetry transformations, SU{3)c x SU{2)l x U{1)y x U{1)x- 

Dim-5 and dim-6 operators with ipip will contribute to thermalization of ip and -0^: 

dim-5: -^^^pH^H , -^^^tpx^ X , ^Ipa^'-^B^, , A.^a'^-^x^, (9.1) 
dim-6: -l^^J^iJJJ^^f , ^V^7^V [^^M - (D^i^tjij] , etc. (9.2) 



^°In case local dark symmetry is spontaneously broken, dark matter candidates may decay via nonrenormalizable 
operators. This will be discussed in detail in a separate publication [9]. 

^^One can refer to Ref. [76, 77] for dim-5 and dim-6 operators in the SM case. 
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where / is the SM chiral fermion field. The first two dim-5 operators above contribution to ^/''i/' 
HH\XX\ whose cross section is estimated as 



1 



in the hmit ^ rriH^rax which is legitimate in our model. This is far less than the cross section 
into a pair of massless dark photon, Eq. (3.4), derived in Sec. 3.1 even if the dark gauge coupling is 
very small. The dim-6 and dim-7 operators will be even smaller. Therefore 

Dim-6 operators including X'^ X will contribute to thermalization of X and 



dim — 6 : 



M2 



X^D^X) - {D^X^)X 



1 



X\D^X) - {D^X^)x] \h^D^H - {D^'H^)H 



etc. 



(9.3) 
(9.4) 



where are the dim-4 gauge invariant SM operators appearing in the SM lagrangian. The cross 
section for annihilation of X and X'^ from dim-6 operators is estimated as 



a{XX'^ SM particles) 



m 



X 



which is totally negligible compared with the annihilation through renormalizable Higgs portal inter- 
action involving Xhx- 

We do not show l/M|-suppressed dim-6 operators for ip decay into XIH, since it is far subdominant 
than the dim-6 operators generated by the virtual RH neutrinos which we already studied in Sec. 5. 



10 Conclusion 

In this paper, we showed that, if the dark matter stability is guaranteed by unbroken local dark 
symmetry with nonzero dark charge, renormalizable portal interactions of the RH neutrinos and SM 
Higgs {H^H) fix the minimal field contents of dark sector (a scalar X and fermion tp as well as mass- 
less dark photon X^) and allow very rich physics without conflicting with various phenomenological, 
astrophysical and cosmological constraints coming from the existence of massless dark photon. 

The unbroken local dark symmetry is very strongly constrained by small and large scale structure 
formation, requiring the dark fine structure constant to be 

ax < 10-^ - 10-' (10.1) 

for O(10^~^)GeV scale mass of dark matter. On the other hand, the dark interaction can be the 
solution to the core/cusp and "too big to fail" problems of small scale structure in collisionless CDM 
scenario. The smallness of dark interaction could cause a danger of dark matter over-abundance. In 
our model, this potential danger is removed by RH neutrino portal which allows ip to decay to X, and 
Higgs portal which allows efficient dilution of the stable dark matter X to get a right amount of dark 
matter relic density. All these nice features are consequences of local dark gauge symmetry, and the 
assumption that all the SM singlet operators being portals to the dark sector. The RH neutrino portal 
also allows production of dark sector asymmetry as leptogenesis in type-I seesaw model does in the 
visible sector, but the dark sector asymmetry eventually disappears as ip decays and does not play any 
significant role. However, one should note that in our scenario eventual relic density of dark matter can 

^^There are no gauge invariant dim-5 operators involving X. 
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be determined by thermal or non-thermal freeze-out, depending on the temperature when ■0 decays to 
X. This allows wider range of dark matter annihilation cross section. Additionally, depending on ax 
and the mass of ip, the decays ip and ip can be the origin of the present baryon number asymmetry 
irrespective of the possible asymmetry between -0 and ip. 

The Higgs portal interaction to dark scalar X cures instability problem of the SM Higgs potential 
by loop effect. So, by introducing large non-minimal gravitational couplings to scalar fields, it becomes 
possible to realize Higgs inflation whose high enough reheating temperature sets the initial condition 
for leptogenesis in our model. The portal interactions also make the dark sector be accessible by 
direct and/or indirect searches, which are consistent with the current bounds from various terrestrial 
experiments and observations in the sky. 

It turned out that the contribution of dark photon to the radiation density at present is about 6% 
of the energy density of a massless neutrino. It is rather small, and still consistent with the present 
observation within 2-a error. The smallness is originated from the fact that dark photon couples only 
to the dark sector fields and dark matter is decoupled from SM particles before QCD-phase transition. 

Our model is a sort of the minimal model which has inflation, leptogenesis, absolutely stable dark 
matter and dark radiation, as well as seesaw mechanism for neutrino masses and mixings, modulo 
some variations described at the end of Sec. 8. The basic principles for the model building were the 
local gauge symmetry working on the dark sector too, and the assumption of the SM singlet operators 
being portals to the dark sector. From these simple principles, we could derive very rich physics 
results that are fully consistent with all the observations made so far. It is interesting to note that the 
Higgs property measurements will strongly constrain our model, since we predict that the Higgs signal 
strength should be equal to "1" for all the decay channels of Higgs boson. Especially the invisible 
Higgs decay into a pair of CDM is kinematically forbidden. 

Our model could be extended in various directions. First of all, one can consider the case where the 
dark gauge symmetry is spontaneously broken. Depending on the local dark symmetry being broken 
fully or partially in case of non-Abelian dark gauge symmetry, one would have interesting phenomeno- 
logical consequences in both dark matter and Higgs sectors. Secondly, it would be straightforward to 
extend our model to super symmetric one. The qualitative features of SUSY models will be similar 
to this letter. However there are two types of CDM. One in the MSSM sector, and the other in the 
other sector. Therefore one could consider the i2-parity in the MSSM sector is broken (either lepton 
number or baryon number). Then the Gx-charged hidden sector dark matter will make the one in the 
universe. These issues will be addressed in the future publication. 

Note Added 

After we submit this paper on the archive, we received a new result by the Planck Collaboration on 
the effective number of neutrino species to be N^^ = 3.30 it 0.27 at 68% CL [78], which is significantly 
lower than other previous results obtained by WMAP-9, SPT and ACT. It is amusing to notice that 
the new Planck data is in perfect agreement with our prediction (A'eff = 3.106) in the model with 
unbroken U{l)x as well as other cases summarized in Table 2 in Sec. 8. 
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Figure 4. Thermal history of the universe in our model 



A Thermally averaged cross sections 

In this Appendix, we collect the thermally averaged cross sections of dark matter pair annihilations. 
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